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^ ' Abstract 
D 

I'T i ■ Cellular networks are usually modeled by placing the base stations on a grid, with mobile users either 

randomly scattered or placed deterministically. These models have been used extensively but suffer from 
being both highly idealized and not very tractable, so complex system-level simulations are used to evaluate 
coverage/outage probability and rate. More tractable models have long been desirable. We develop new general 
models for the multi-cell signal-to-interference-plus-noise ratio (SINR) using stochastic geometry. Under very 
general assumptions, the resulting expressions for the downlink SINR CCDF (equivalent to the coverage 
probability) involve quickly computable integrals, and in some practical special cases can be simplified to 
common integrals (e.g., the Q-function) or even to simple closed-form expressions. We also derive the mean 
rate, and then the coverage gain (and mean rate loss) from static frequency reuse. We compare our coverage 
predictions to the grid model and an actual base station deployment, and observe that the proposed model 
is pessimistic (a lower bound on coverage) whereas the grid model is optimistic, and that both are about 

o 

equally accurate. In addition to being more tractable, the proposed model may better capture the increasingly 
opportunistic and dense placement of base stations in future networks. 

X 

c5 ■ I. Introduction 

Cellular systems are now nearly universally deployed and are under ever-increasing pressure to 
increase the volume of data they can deliver to consumers. Despite decades of research, tractable 
models that accurately model other-cell interference (OCI) are still unavailable, which is fairly re- 
markable given the size of the industry. This deficiency has impeded the development of techniques 
to combat OCI, which is the most important obstacle to higher spectral efficiency in today's cellular 
networks, particularly the dense ones in urban areas that are under the most strain. In this paper we 
develop accurate and tractable models for downlink capacity and coverage, considering full network 
interference. 
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A. Common Approaches and Their Limitations 

A tractable but overly simple downlink model commonly used by information theorists is the Wyner 
model HI, 0, J3]|, which is typically one-dimensional and presumes a unit gain from each base station 
to the tagged user and an equal gain that is less than one to the two users in the two neighboring 
cells. This is a highly inaccurate model unless there is a very large amount of interference averaging 
over space, such as in the uplink of heavily-loaded CDMA systems [4J. This philosophical approach 
of distilling other-cell interference to a fixed value has also been advocated for CDMA in [[51 and 
used in the landmark paper ||6), where other-cell interference was modeled as a constant factor of the 
total interference. For cellular systems using orthogonal multiple access techniques such as in LTE 
and WiMAX, the Wyner model and related mean-value approaches are particularly inaccurate, since 
the SINR values over a cell vary dramatically. Nevertheless, it has been commonly used even up to the 
present to evaluate the "capacity" of multicell systems under various types of multicell cooperation 
0, (HI, [0. Another common analysis approach is to consider a single interfering cell ifTOl . In the two 
cell case, at least the SINR varies depending on the user position and possibly fading, but naturally 
such an approach still neglects most sources of interference in the network and is highly idealized. A 
recent discussion of such models for the purposes of base station cooperation to reduce interference is 
given in [fTTIl . That such simplified approaches to other-cell interference modeling are still considered 
state-of-the-art for analysis speaks to the difficulty in finding more realistic tractable approaches. 

On the other hand, practicing systems engineers and researchers in need of more realistic models 
typically model a 2-D network of base stations on a regular hexagonal lattice, or slightly more simply, 
a square lattice, see Fig. [2] Tractable analysis can sometimes be achieved for a fixed user with a small 
number of interfering base stations, for example by considering the "worst-case" user location - the 
cell corner - and finding the signal-to-interference-plus-noise ratio (SINR) lfT2l . [fT3l . The resulting 
SINR is still a random variable in the case of shadowing and/or fading from which performance 
metrics like (worst-case) average rate and (worst-case) outage probability relative to some target rate 
can be determined. Naturally, such an approach gives very pessimistic results that do not provide much 
guidance to the performance of most users in the system. More commonly, Monte Carlo integrations 
are done by computer, e.g. in the landmark capacity paper (6). Tractable expressions for the SINR are 
unavailable in general for a random user location in the cell and so more general results that provide 
guidance into typical SINR or the probability of outage/coverage over the entire cell must be arrived 
at by complex time-consuming simulations. In addition to being onerous to construct and run, such 
private simulations additionally suffer from issues regarding repeatability and transparency, and they 
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seldom inspire "optimal" or creative new algorithms or designs. It is also important to realize that 
although widely accepted, grid-based models are themselves highly idealized and may be increasingly 
inaccurate for the heterogeneous and ad hoc deployments common in urban and suburban areas, 
where cell radii vary considerably due to differences in transmission power, tower height, and user 
density. For example, picocells are often inserted into an existing cellular network in the vicinity of 
high-traffic areas, and short-range femtocells may be scattered in a haphazard manner throughout a 
centrally planned cellular network. 

B. Our Approach and Contributions 

Perhaps counter-intuitively, this paper addresses these long-standing problems by introducing an 
additional source of randomness: the positions of the base stations. Instead of assuming they are 
placed deterministically on a regular grid, we model their location as a homogeneous Poisson point 
process of density A. Such an approach for BS modelling has been considered as early as 1997 |[T4l . 
031, lfT6l but the key metrics of coverage (SINR distribution) and rate have not been determined!! The 
main advantage of this approach is that the base station positions are all independent which allows 
substantial tools to be brought to bear from stochastic geometry; see [fTSl for a recent survey that 
discusses additional related work, in particular lfT9"l . Il20l . [|2T|. Although BS's are not independently 
placed in practice, the results given here can in principle be generalized to point processes that model 
repulsion or minimum distance, such as determinantal and Matern processes [|22l . Il23l . The mobile 
users are scattered about the plane according to some independent homogeneous point process with a 
different density, and they communicate with the nearest base station while all other base stations act 
as interferers, as shown in Fig. [TJ 

From such a model, we achieve the following theoretical contributions. First, we are able to derive 
a general expression for the probability of coverage in a cellular network where the interference 
fading/shadowing follows an arbitrary distribution. The coverage probability is the probability that a 
typical mobile user is able to achieve some threshold SINR, i.e. it is the complementary cumulative 
distribution function (CCDF). This expression is not closed-form but also does not require Monte 
Carlo methods. The coverage is then derived for a number of special cases, namely combinations of 
(i) exponentially distributed interference power, i.e. Rayleigh fading, (ii) path loss exponent of 4, and 
(iii) interference-limited networks, i.e. thermal noise is ignored. These special cases have increasing 
tractability and in the case that all three simplifications are taken, we derive a remarkably simple 



'The paper 1171 was made public after submission of this paper and contains some similar aspects to the approach in this paper. 



4 



formula for coverage probability that depends only on the threshold SINR. We compare these novel 
theoretical results with both traditional (and computationally intensive) grid-based simulations and 
with actual base station locations from a current cellular deployment in a major urban area. We see 
that over a broad range of parameter and modeling choices our results provide a reliable lower bound 
to reality whereas the grid model provides an upper bound that is about equally loose. In other words, 
our approach, even in the case of simplifying assumptions (i)-(iii), appears to not only provide simple 
and tractable predictions but also accurate ones. 

Next, we derive the mean achievable rate in our proposed cellular model under similar levels of 
generality and tractability. The two competing objectives of coverage and rate are then explored 



analytically through the consideration of frequency reuse, which is used in some form in nearly 
all cellular systems^ to increase the coverage or equivalently the cell edge rates. Our expressions for 
coverage and rate are easily modified to include frequency reuse and we find the amount of frequency 
reuse required to reach a specified coverage probability, as well as seeing how frequency reuse degrades 
mean rate by using the total bandwidth less efficiently. 

II. Downlink System Model 

The cellular network model consists of base stations (BSs) arranged according to some homogeneous 
Poisson point process (PPP) $ of intensity A in the Euclidean plane. Consider an independent collection 
of mobile users, located according to some independent stationary point process. We assume each 
mobile user is associated with the closest base station; namely the users in the Voronoi cell of a BS 
are associated with it, resulting in coverage areas that comprise a Voronoi tessellation on the plane, as 
shown in Fig. [Q A traditional grid model is shown in Fig. [2] and an actual base station deployment in 
Fig. [3l The main weakness of the Poisson model is that because of the independence of the PPP, BSs 
will in some cases be located very close together but with a significant coverage area. This weakness is 
balanced by two strengths: the natural inclusion of different cell sizes and shapes and the lack of edge 
effects, i.e. the network extends indefinitely in all directions. The models are quantitatively compared 
in Section |V] 

The standard power loss propagation model is used with path loss exponent a > 2. As far as random 
channel effects such as fading and shadowing, we assume that the tagged base station and tagged user 
experience only Rayleigh fading with mean 1, and constant transmit power of Then the received 



2 Even cellular systems such as modern GSM and CDMA networks that claim to deploy universal frequency reuse still thin the 
interference in time or by using additional frequency bands - which is mathematically equivalent to thinning in frequency. 
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power at a typical node a distance r from its base station is hr~ a where the random variable h 
follows an exponential distribution with mean If/j,, which we denote as h ~ exp(/x). Note that other 
distributions for h can be considered using Prop. 2.2 of Il24ll but with some loss of tractability. The 
interference power follows a general statistical distribution g that could include fading, shadowing, 
and any other desired random effects. Simpler expressions result when g is also exponential and these 
are given as special cases. Lognormal interference is considered numerically^, and we see although 
it degrades coverage it does not significantly affect the accuracy of our analysis. Because of these 
random channel effects, in our model not all users will be connected to the base station capable of 
providing the highest SINR. All results are for a single transmit and single receive antenna, although 
future extensions to multiple such antennas are clearly desirable. 

The interference power at the typical receiver I r . is the sum of the received powers from all other 
base stations other than the home base station and is treated as noise in the present work. There is no 
same-cell interference, for example due to orthogonal multiple access within a cell. The noise power 
is assumed to be additive and constant with value a 2 but no specific distribution is assumed in general. 
The snr = ^2 is defined to be the received snr at a distance of r = 1. All analysis is for a typical 
mobile node which is permissible in a homogeneous PPP by Slivnyak's theorem ll23l . 



III. Coverage 

This is the main technical section of the paper, in which we derive the probability of coverage in 
a downlink cellular network at decreasing levels of generality. The coverage probability is defined as 

p e (T,\,a) ^P[sinr>T], (1) 

and can be thought of equivalently as (i) the probability that a randomly chosen user can achieve a 
target SINR T, (ii) the average fraction of users who at any time achieve SINR T, or (iii) the average 
fraction of the network area that is in "coverage" at any time. The probability of coverage is also 
exactly the CCDF of SINR over the entire network, since the CDF gives P[sinr < T]. 

Without any loss of generality we assume that the mobile user under consideration is located at the 
origin. A user is in coverage when its SINR from its nearest BS is larger than some threshold T and 
it is dropped from the network for SINR below T. The SINR of the mobile user at a random distance 

r from its associated base station can be expressed as 

hr~ a 

SINR = , (2) 

cH + I r 

3 Shadowing is neglected between the tagged BS and user since it can fairly easily be overcome with even slow power control. In 
this case the transmit power would be simply l/g/-t and treated as a constant over the shadowing time-scale. 
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where 

I r = 9iRj a (3) 

is the cumulative interference from all the other base stations (except the tagged base station for the 
mobile user at o denoted by b Q ) which are a distance Ri from the typical user and have fading value 

9i- 

A. Distance to Nearest Base Station 

An important quantity is the distance r separating a typical user from its tagged base station. Since 
each user communicates with the closest base station, no other base station can be closer than r. In 
other words, all interfering base stations must be farther than r. The probability density function (pdf) 
of r can be derived using the simple fact that the null probability of a 2-D Poisson process in an area 
A is exp(-XA). 

P[r > R] = P[No BS closer than R] (4) 
= e- x * R \ (5) 



Therefore, the cdf is P[r < R] — F r (R) = 1 — e XnR2 and the pdf can be found as 

Mr) - ^ (6) 
dr 

= e- Xnr2 27T\r. (7) 

Meanwhile, the interference is a standard M/M shot noise [|25l . Il22ll . E6l created by a Poisson 
point process of intensity A outside a disc at center o and of radius r, for which some useful results 
are known and applied in the sequel. 

B. General Case and Main Result 

We now state our main and most general result for coverage probability from which all other results 
in this section follow. 

Theorem 1: The probability of coverage of a typical randomly located mobile user in the general 
cellular network model of Section HI] is 



POO 

Pc (T, A, a) = VTA / e -*At0(T,a)-MT^W* dVj 

Jo 



(8) 



where 

1(nT\l r , i 

(9) 



2(uTU 

P(T, a) = — ^— ^— E gs (r(-2/a, jiTg) - T(-2/a)) 



a 
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and the expectation is with respect to the interferer's channel distribution g. Also, F(a, x) = t a ~ 1 e~ t dt 
denotes the incomplete gamma function, and T(x) = J °° t x ~ l e~ l dt the standard gamma function. 

Proof: Conditioning on the nearest BS being at a distance r from the typical user, the probability 
of coverage averaged over the plane is 

p c (T, A, a) = E r [P[sinr > T \ r]} 

= / P[sinr > T | r]f r (r)dr 

Jr>0 



(a) 



P 



r>0 



r>0 



hr 



a 2 + I r 



> T 



'F[h> Tr a (a 2 + I r ) \ r}2n\rdr. 



The distribution f r (r) and hence (a) follows from Subsection IIII-Al Using the fact that h ~ exp(fi), 
the coverage probability can be expressed as 

F[h > Tr a (a 2 + 7 r ) | r] = E /r [¥[h > Tr a (a 2 + I r ) \ r, I r 

= E /r [exp(-fiTr a (a 2 + I r )) \ r] 



(10) 



where £i r (s) is the Laplace transform of random variable I r evaluated at s conditioned on the distance 
to the closest BS from the origin. This gives a coverage expression 



p c (T, A, a) 



C Ir (fiTr a )27r\rdr. 



(ID 



r>0 



Defining Ri as the distance from the ith interfering base station to the tagged receiver and g,i as the 
interference channel coefficient of arbitrary but identical distribution for all i, using the definition of 
the Laplace transform yields 

£ /r («) =E /p [e- / -]=E. lft [eK P (-« ^R- a )\ 

iS$\{6 } 



E 



n ex p(- s ^^ : 

iS*\{6 } 



(a) 



J] E g [exp(-sgR7")} 
■ie<s>\{b } 



exp [exp(— sgv a )]) vdv J , 



(12) 



s 



where (a) follows from the i.i.d. distribution of Qi and its further independence from the point process 
$, and the last step follows from the probability generating functional (PGFL) G31 of the PPP, which 
states for some function f(x) that E [n^e* f( x )] = ex P (~ ^ IrsO- ~ f( x ))^ x ) ■ The integration limits 
are from r to oo since the closest interferer is at least at a distance r. Let f(g) denote the PDF of g. 
Plugging in s = fxTr a , and swapping the integration order gives, 

C Ir (^Tr a ) = exp (-2vrA £° (T°(l - e^^vdi^j f(g)dg) . 

The inside integral can be evaluated by using the change of variables v~ a — > y, and the Laplace 
transform is 

£ /r (/iTr«) = exp (w _ 2 ^J)^ 2 ^ [T ^ 2 /a^Tg) - T(-2/a)] f(g)dg 

Combining with (TTTT) . and using the substitution r 2 — > v, we obtain the result. ■ 
In short, Theorem 1 gives a general result for the probability of achieving a target SINR T in the 
network. It is not closed-form but the integrals are fairly easy to evaluate. We now turn our attention 
to a few relevant special cases where significant simplification is possible. 

C. Special Cases: Interference Experiences General Fading 

The main simplifications we will now consider in various combinations are (i) allowing the path 
loss exponent a = 4, (ii) an interference-limited network, i.e. 1/cr 2 — > oo, which we term "no noise" 
and (iii) interference fading power g ~ exp(/i) rather than following an arbitrary distribution! In this 
subsection we continue assume the interference power follows a general distribution, so we consider 
two special cases corresponding to (i) and (ii) above. 

1 ) General Fading, Noise, a = 4: First, if a = 4, Theorem 1 admits a form that can be evaluated 
according to 

oo r~ / 2 

-ax-bx 2 - I" I a 



e-- e -«d* = ,/ r x P (-j)QI 71 ), (13) 



where Q(x) = ^= exp(— y 2 /2)dy is the standard Gaussian tail probability. Setting a = tt\(3(T, a) 
and b = fiTa 2 = T/snr gives 

*(T,A,4) = ^expf^^W^MY (14) 
1 ^ V^T^ V 4T/snr J H \ ^XTJ^k J 

Therefore, given the numerical calculation of (3(T, 4) for a chosen interference distribution, the 
coverage probability can be found in quasi-closed form since Q(x) can be evaluated nearly as easily 
as a basic trigonometric function by modern calculators and software programs. 



4 The interference power is also attenuated by the path loss so the mean interference power for each base station is less than the mean 
desired power, by definition, even though the fading distributions have the same mean y,, which is a proxy for the transmit power. 
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2) General Fading, No Noise, a > 2: In most modern cellular networks thermal noise is not an 
important consideration. It can be neglected in the cell interior because it is very small compared 
to the desired signal power (high SNR), and also at the cell edge because the interference power is 
typically so much larger (high INR). If a 2 — > (or transmit power is increased sufficiently), then 
using Theorem 1 it is easy to see that 



In the next subsection, we show that (fl~4"l) does in fact reduce to (fl~5T) as a 2 — > 0, which is not obvious 
by inspection. 

It is interesting to note that in this case the probability of coverage does not depend on the base station 
density A. It follows that both very dense and very sparse networks have a positive probability of cov- 
erage when noise is negligible. Intuitively, this means that increasing the number of base stations does 
not affect the coverage probability, because the increase in signal power is exactly counter-balanced by 
the increase in interference power. This matches empirical observations in interference-limited urban 
networks as well as predictions of traditional, less-tractable models. In interference-limited networks, 
increasing coverage probability typically requires interference management techniques, for example 
frequency reuse, and not just the deployment of more base stations. Note that deploying more base 
stations does allow more users to be simultaneously covered in a given area, both in practice and 
under our model, because we assume one active user per cell. 

3) General Fading, Small but Non-zero Noise: A potentially useful low noise approximation of 
the success probability can be obtained that is more easily computable than the constant noise power 
expression and more accurate than the no noise approximation for a 2 ^ 0. Using the expansion 
exp(— x) = 1 — x + o{x), x — > it can be found after an integration by parts of © that 



For the special case of a = 4, it is not immediately obvious that (fl~4)) is equivalent to (fl~5l) as a 2 — > 0, 
but indeed it is true. It is possible to write (fT4l) as 



p c (T, A, a) 



1 



P(T,aY 



(15) 




(16) 




(17) 



where x 



V2b 



a 



and a, b as before. The series expansion of Q(x) for large x is 




(18) 



which means that 




(19) 
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which allows simplification of ( TT71) to (fl~5l) for the case of no noise. 

D. Special Cases: Interference is Rayleigh Fading 

Significant simplification is possible when the interference power follows an exponential distribution, 
i.e. interference experiences Rayleigh fading and shadowing is neglected. We give the coverage 
probability for this case as Theorem 2. 

Theorem 2: The probability of coverage of a typical randomly located mobile user experiencing 
exponential interference is 

poo 

Pc (T, A, a) = n\ e -*Mi+ P (T.«))-/^W> dUj (20) 
Jo 

where 

poo 1 

p(T,a)=T 2 M —— -^—du. (21) 

J T -2/ a 1 + U a ' 2 

Proof: The proof is a special case of Theorem [H but however lends to much simplification. The 
proof is provided in Appendix [A] ■ 
We now consider the special cases of no noise and a = 4. 

1 ) Exponential Fading, Noise, a = A: When a = 4, using the same approach as in (fl"3l) . we get 

vrlA [(\txk(T)) 2 \ „ ( XttkCT) \ 



where «(T) = 1 + p(T,4) = 1 + Vf(n/2 - arctan(l/v / T)). 

This expression is quite simple and is practically closed-form, requiring only the computation of a 
simple Q(x) value. 

2) Exponential Fading, No Noise, a > 2: In the no noise case the result is very similar to general 
fading in appearance, i.e. 

Pc(T, A, a) = — (23) 
1 + p(T, a) 

with p(T, a) being faster and easier to compute than the more general expression (3(T, a). When the 
path loss exponent a = 4, the no noise coverage probability can be further simplified to 

p c (T, A, 4) = = — . (24) 

1 + VT(tt/2 - arctan(l/v / T)) 

This is a remarkably simple expression for coverage probability that depends only on the SIR threshold 
T, and as expected it goes to 1 for T — > and to for T — > oo. For example, if T = 1 (0 dB, which 
would allow a maximum rate of 1 bps/Hz), the probability of coverage in this fully loaded network 
is 4(4 + 7r) _1 = 0.56. This will be compared in more detail to classical models in Section |V] A small 
noise approximation can be performed identically to the procedure of Section ITII-C3I with 1 + p(T, a) 
replacing f3(T,a) in (fT6l) . 
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IV. Average Achievable Rate 

In this section, we turn our attention to the mean data rate achievable over a cell. Specifically we 
compute the mean rate in units of nats/Hz (1 bit = ln(2) = 0.693 nats) for a typical user where 
adaptive modulation/coding is used so each user can set their rate such that they achieve Shannon 
bound for their instantaneous SINR, i.e. ln(l + sinr). Interference is treated as noise which means 
the true channel capacity is not achieved, which would require a multiuser receiver 11271 . 11281 . [|29l , 
but future work could relax this constraint within the random network framework, see e.g. [|30l . OTTl . 
In general, almost any type of modulation, coding, and receiver structure can be easily treated by 
adding a gap approximation to the rate expression, i.e. r — > ln(l + sinr/G) where G > 1 is the gap. 
The technical tools and organization are similar to Section [HI] so the discussion will be more concise. 
The results are all for exponentially distributed interference power but general distributions could be 
handled as well following the approach of Theorem \T\ and techniques from ll24l . 



A. General Case and Main Result 

We begin by stating the main rate theorem that gives the ergodic capacity of a typical mobile user 
in the downlink. 

Theorem 3: The average ergodic rate of a typical mobile user and its associated base station in the 
downlink is 

r(A,a) = E[ln(l + sinr)] (25) 
= f e~ nXr2 [ e-° Va < ei - 1 >£j r (/ir a (e* - l))dt27rArdr, (26) 

Jr>0 Jt>0 

where 

£ /r (/irV - 1)) = exp f-vrArV - 1) 2 / Q T T^T**^) ' (27) 

Proof: The proof is provided in Appendix El ■ 
The computation of r in general requires three numerical integrations. 

B. Special Case: a = 4 

For a = 4 the mean rate simplifies to 

r (A,4)= [ / e - CT2 ^ 4 ( et - 1 )e- wAr H 1+ ^( w/2 " arctan(1/v/?3T) ))27rArdrdt. 



<>0r>0 

e -a 2 ^ 2 (e t -l)/M) 2 e -t>(l+v / ?=T(7r/2-arctan(l/v / i^T))) dt;dt 




t>0v>0 
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Using dH), 



t>o 



b(t) ^ \m J ^ \VW) J 



^( A > 4 ) = / Jl7IT ex P ( 7I77T ) Q I "9= ) d ^ ( 28 ) 



where a(t) = 1 + (tt/2 - arctan(l/ v / e F ^T)) and 6(f) = a 2 /i(e* - 1)/(ttA) 2 . The final 

expression (1281) be evaluated numerically with one numerical integration (presuming an available look 
up table for Q(x)). 

C. Special Case: No Noise 

When noise is neglected, a 2 — > 0, so 

r(X,a) = [ [ exp f -ix\r 2 ( 1 + (e* - l) 2/a f ^dx ) ) 2vrArdrdt 

Jr>0 Jt>0 V V J {e t-i)-y« 1 + x a ' 2 J J 

exp l-v (l + (e* - l) 2/a / ^jdx J J dwdt 



t>0 ir>0 V V J(e t -l)-2/a 1 + 

i>0 1 + \ e L ) 1 J(e*-l)-2/<« 1+W 2 

a quantity that again does not depend on A. As in the case of coverage, increasing the base station 
density does not increase the interference-limited ergodic capacity per user in the downlink because 
the distance from the mobile user to the nearest base station and the average distance to the nearest 
interferer both scale like 0(A -1 / 2 ), which cancel. Note, however, that the overall sum throughput and 
area spectral efficiency of the network do increase linearly with the number of base stations since 
the number of active users per area achieving rate r is exactly A, assuming that the user density is 
sufficiently large such that there is at least one mobile user per cell. 
In the particular case of a = 4 in conjunction with no noise, 



e 



oo 

t _ i\2/« 



r°° i / 

/ rrdx = Ve* — 1 ( 7r/2 — arctanfl/vV — 1 

y fe *_n- 2 /a 1 + a-/ 2 V 1 y 1 



'( e *_l)-2/c 

so the mean rate is expressed to a single simple numerical integration that yields a precise scalar 

r(A,4)=/ - dt re 1.49nats/sec/Hz. (30) 

Jt>o 1 + Ve^T (n/2 - arctan(l/v / i F ^T)) 

In other words, our model predicts that the no noise limit for mean downlink rate in a cellular system 
with Rayleigh fading is 2.15 bps/Hz if a = 4. 
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V. Validation of the Proposed Model 

Now that we have developed expressions for the coverage and mean rate for cellular networks, it is 
important to see how these analytical results compare with the widely accepted grid model. Further, 
we were able to obtain precise coordinates for base stations over a large urban area from a major 
service provider, and we compare our results to the coverage predicted by those locations as well 
(which are neither a perfect grid nor Poisson). Intuitively, we would expect the Poisson model to 
give pessimistic results compared to a planned deployment due to the strong interference generated 
by nearby base stations. The grid model is clearly an upper bound since a perfectly regular geometry 
is in fact optimal from a coverage point of view lfT6ll . An additional source of optimism in the grid 
model is the customary neglect of background interference from outer tier base stations. We see in 
Section IV-AI that the latter effect is not very significant, however. 



A. The Grid Model and An Actual BS Deployment 

A periodic grid is typically used in prior work to model the base station locations. We use a square 
lattice for notational simplicity but a hexagonal one can also be used all results will only differ by 
a very small constant. We consider a home base station located at the origin and N interfering base 
stations located in square tiers around the home base station. Each tier is a distance 2R from the 
previous tier, i.e. each base station coverage area is a 2R x 2R square, and so any user within a 
distance R of a base station is guaranteed to be covered by it. The base station density in this case is 
1/AR 2 base stations per unit area. A two tier example with N = 24 is shown in Figure |2] The SINR 
for a regular base station deployment becomes 

hr~ a 

SINR = — u —, (31) 

In + O" 2 

where r u = \/ x 2 u + y 2 with x u ~ U[—R,R] and y u ~ U[—R,R]. The channel fading power is still 
h ~ exp(/i) as in previous sections. The interference to the tagged user is now 

N 

4 = 5>rT Q (32) 

i=l 

where r\ = a/ (x{ — x u ) 2 + — y u ) 2 is the distance seen from interfering base station % and gi its 
observed fading power. The probability of coverage is 

p c (T, a) = P[sinr >T]= F[h > r a u T{I u + a 2 )], (33) 

which is no different in principle than ©, but due to the structure of I u it is difficult to proceed 
analytically, and so numerical integration is used to compare with the results of Sections [III] and 
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ITVl One important difference between the behavior of grid and random BS models are the allowed 
extremes on the distances of the tagged and interfering base stations. In a grid model, there is always a 
base station within a specified distance R and never an interfering one closer than R. In the proposed 
model, two base stations can be arbitrarily close together and hence there is no lower bound on R, so 
both the tagged and an interfering base station can be arbitrarily close to the tagged user. The only 
constraint is that the interfering base station must be farther than the tagged one, or else a handoff 
would occur. 

We have also obtained the coordinates of a current base station deployment by a major service 
provider in a relatively flat uniform urban area. This deployment stretches over an approximately 
100 x 100 km square, and we show a zoom of the middle 40 x 40 km in Fig. [3l In this figure, the cell 
boundaries correspond to a Voronoi tessellation and hence are only a function of Euclidean distance, 
whereas in practice other factors might determine the cell boundaries. Clearly this is only a single 
deployment and further validation should be done. However, we strongly suspect that deployments in 
many cities follow an even less regular topology due to irregular terrain such as large hills and water 
features and/or high concentrated population centers. It seems such scenarios might be even better 
suited to a random spatial model that the example provided here. 

B. Coverage Comparison 

In Fig. HI we compare the traditional square grid model to the random PPP base station model. The 
plot gives the probability that a given sinr target T on the x-axis can be achieved, i.e. it gives the 
complementary cumulative distribution function (CCDF) of sinr, i.e. P[sinr > T). Both N = 8 and 
N = 24 are used, and it can be seen that the N = 8 case is only slightly more optimistic as opposed 
to N = 24, at least for a = 4 (the gap increases slightly for smaller a). The curves all exhibit the 
same basic shape and as one would expect, a regular grid provides a higher coverage area over all 
possible sinr targets. A small (< 1 dB) gap is seen between the snr = 10 and snr — > oo cases, which 
confirms that noise is not a very important consideration in dense cellular networks, which are known 
to be interference-limited. Therefore, we neglect noise in the ensuing plots. 

In Fig. [5] we now compare the three different base station location models with exponential (Rayleigh 
fading) interference. The random BS model is indeed a lower bound and the grid model an upper 
bound. The random BS model appears no worse than the grid model in terms of accuracy and may be 
preferable from the standpoint that it provides conservative predictions, as well as being much more 
analytically tractable. The Poisson BS model becomes more accurate at lower path loss exponents. 
There are two reasons for this. First, the PPP models distant interference whereas a 1 or 2 tier grid 
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model does not; and the the interference of far-off base stations is more significant for small a. Second, 
since a weakness of the Poisson model is the artificially high probability of a nearby and dominant 
interfering base station, at lower path loss exponents, perhaps counter-intuitively, such an effect is less 
corrupting because a dominant base station contributes a lower fraction of the total interference due 
to the slower attenuation of non-dominant interferers. 

Next we consider the effect on lognormal interference, which is common in cellular networks. 
Whereas shadowing to the desired base station can be overcome with power control (or macrodiversity, 
not considered here) the interference remains lognormal. We assume the shadowing is given by a value 
10^> where X ~ N(£, k 2 ) and £ and k are now in dB. We normalize £ to be the same as for the 
exponential case and consider various values of k in Figs. |6] and [7J Fig. [6] shows the extent to which 
lognormal interference increases the coverage probability in our model, whereas Fig. [7J shows that our 
model still reasonably tracks a real deployment. It may seem counterintuitive that increasing lognormal 
interference increases the coverage probability, the reason being that cell edge users have poor mean 
SINR (often below T), and so increasing randomness gives them an increasing chance of being in 
coverage. It also implies that SINR-aware scheduling, which is not considered here, might be able to 
significantly increase coverage. 

VI. Frequency Reuse: Coverage vs. Rate 

Cellular network operators must provide at least some coverage to their customers with very high 
probability. For example, sinr = 1 might be a minimal level of quality needed to provide a voice call. 
In this case, for a = 4 we can see from Fig. @] that the grid model gives a success probability of 
about 0.7 and the PPP model predicts 0.53. Clearly, neither is sufficient for a commercial network, 
so cellular designers must find a way to increase the coverage probability. Assuming the network is 
indeed interference-limited, a common way to do this is to reduce the number of interfering base 
stations. This can be done statically through a planned and fixed frequency reuse pattern and/or cell 
sectoring, or more adaptively via a reduced duty cycle in time (as in GSM or CDMA voice traffic), 
fractional frequency reuse, dynamic bandwidth allocation, or other related approaches 11321 . |[33l . More 
sophisticated interference cancellation/suppression approaches can also be used, potentially utilizing 
multiple antennas. In this paper, we restrict our attention to straightforward per-cell frequency reuse. 

In frequency reuse, the reuse factor 6 > 1 determines the number of different frequency bands used 
by the network, where just one band is used per cell. For example, if 5 = 2 then the square grid of 
Fig. [2] can assign the top row base stations frequencies 1, 2, 1, 2, 1, and then the second row 2, 1, 
2, 1,2, and so on. In this way interfering base stations are now separated by a distance 2^/2R rather 
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than 2R. Larger values of S monotonically decrease the interference, e.g. 5 = 4 allows a base station 
separation of AR in the square grid model. 

The PPP BS model also allows for interference thinning, but instead of a fixed pattern (which is not 
possible in a random deployment) we assume that each base station picks one of 5 bands at random. 
A visual example is given in Fig. [8] for 5 = 4. The model suffers from the fact that adjacent base 
stations may simultaneously use the same frequency even for large 5. However, it is not clear this is 
any worse of a model than the idealistic grid model with pre-planned frequency reuse, especially for 
4G OFDMA-based networks that will use dynamic frequency allocation with very aggressive overall 
frequency reuse, so that each subcarrier may appear from a birds-eye view of the network to be 
allocated almost randomly at any snapshot in time. 



A. Increasing Coverage via Frequency Reuse 

First, we consider the effect of random frequency reuse on the coverage probability. 
Theorem 4: If 5 frequency bands are randomly allocated to the cells, then the coverage probability 
with exponentially distributed interference power is equal to 

p e {T, A, a, 5)=tt\J^ exp f-nXv (l + ^p{T, a)^j - fiTa 2 v a/2 j dv. (34) 

Proof: A typical mobile user at the origin o would be served by its closest BS from the complete 
point process $. Call this distance r, then as in Section IIII-AI r is Rayleigh distributed with PDF 
f r (r) = \2irr exp(—\nr 2 ) with the interferers located outside r. The interfering BSs which transmit 
in the same frequency band are a thinned version of the original PPP and have a density X/5. Since 
a thinned version of a PPP is again a PPP, the rest of the proof exactly follows Theorem [2l ■ 

Observe that as the number of frequency bands 5 — > oo a coverage limit is reached that depends 
only on the noise power. In Fig.[9l the coverage probability is plotted with respect to the sinr threshold 
T for 5 = 2 and 5 = 4 for each of the three BS placement models. 

A cellular operator often wishes to guarantee a certain probability of coverage to its customers. 
For a given blocking/outage probability e, the following corollary for the no noise case provides the 
number of frequency bands that are required. 

Corollary 1: The minimum number of frequency bands needed to ensure an outage probability no 
greater than e is 

5 = K' a)(1 ~ £) l . (35) 
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Proof: For the case of no noise, a 2 = and hence C w {(jTr a ) = 1 which simplifies the coverage 
probability to 

p e (T, A, a, 5) = 3_y, (36) 

from which the result follows by setting this quantity to be equal to 1 — e and requiring it to be an 
integer. ■ 
We now compare these analytical results with central frequency planning (optimal) for the grid model 
and a heuristic approach for the actual base stations. Frequency planning in actual cellular networks 
is a complex optimization problem that depends on the specific geography and data loads and is often 
performed heuristically by the service provider. For nominally reuse 1 systems like LTE, frequency 
reuse can also be done adaptively on a per subcarrier basis but here we consider only static frequency 
reuse. To provide a reasonably fair comparison, we have used a simple centralized greedy frequency 
allocation algorithm that maximizes the distance between cells sharing the same frequency band in the 
actual BS network. Although not "optimal", it provides a reasonable benchmark to compare against 
the performance predictions of PPP model which uses a random frequency allocation. Examples of 
both allocations are shown visually in Fig. [8l Under a random frequency allocation, adjacent cells may 
transmit in the same band with a higher probability as compared to a planned allocation. This leads to 
an increasing gap between the coverage predicted by the PPP model and that of the actual BSs as S 
increases, as per Fig. [9l As expected, central frequency planning also outperforms random frequency 
allocation. This indicates that the proposed model in its current form is not a faithful predictor of 
coverage for large frequency reuse factors. However, a recent work ll34ll uses this approach to study 
related and more general interference management techniques such as fractional frequency reuse. 

B. Frequency Reuse 's Effect on Rate 

The desirable increase in coverage with increasing 5 has to be balanced against the fact that each 
cell can only use l/5th of the available frequencies. In this section we will show that the optimal 5 
from a mean rate point of view is in fact 5=1, i.e. any increase in coverage from frequency reuse is 
paid for by decreasing the overall sum rate in the network. The following general result can be given 
for average rate with frequency reuse, the key observation being that since the bandwidth per cell was 
previously normalized to 1 Hz, now it is i Hz. We assume the SNR per band is unchanged. 

Theorem 5: If 5 frequency bands are randomly allocated to the cells, the average rate of a typical 
mobile user in a downlink is 

T(\,a,5) = - : [ e~ nXr2 f e- CT2 ^ (et - 1) /: /r (/ir Q (e'-l))dr;27rArdr, (37) 

" Jr>0 Jt>0 
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where 

/ Avrr 2 fe* - l) 2/a r°° 1 \ 

br(l*r»V - 1)) = exp - ( % l) / TT~~^/2 d9 ) ■ ( 38 ) 

Proof: The average rate of a typical mobile user is iE[]n(l + sinr)], and the proof proceeds in a 
similar manner to Theorem 0] and Theorem [3] and so is omitted. ■ 
From the above Theorem, the average rate without noise is given by 

r(X,a,5)=[ , f L 1 — j-dt, 

Jt>o $ + (e* - l) 2/a J^t.^-a/a TT^ dx 

and is obviously maximized for 5 = 1. This can also be seen visually in Fig. [TO], which shows the 
average rate as a function of 5 for two different path loss exponents. The next corollary generalizes 
this observation to the case of non-zero noise. 

Corollary 2: The average rate of a typical mobile user r(A, a, S) is maximized for S = 1. 
Proof: Using the substitution r 2 — > 5y in (1371 ), we observe that the integrand decreases with 5, 
hence verifying the claim. ■ 
As in (|3Q|) . the average rate simplifies further for the case of a = 4 and no noise and 

T (XA,S)= / , r dt. 

Jt>o 5 + v 7 ?— I (tt/2 - arctan(l/ v / e F ^T)) 

By using numerical integration, r(A,4, 1) ~ 1.49 nats/sec/Hz, r(A,4, 2) w 1.1 nats/sec/Hz, and 
r(A,4,3) w 0.87 nats/sec/Hz. 

VII. Conclusions 

This paper has presented a new framework for downlink cellular network analysis. It is significantly 
more tractable than the traditional grid-based models, and appears to track (and lower bound) a real 
deployment about as accurately as the traditional grid model (which upper bounds). A final verdict 
on its accuracy will require extensive comparison with further real base station deployments. In view 
of current trends whereby base stations are deployed somewhat opportunistically with ever-increasing 
density, and having variable cell radii, the proposed model may actually become increasingly accurate 
as well as much more tractable. 

Given the number of problems of contemporary interest that require modeling neighboring base 
stations, the possibilities for future work using this model are extensive. An extension to the uplink 
would be desirable. Further extensions to this approach could include random spatial placements of base 
stations that model repulsion, or heterogeneous networks that have both macro and micro/pico/femto 
cells with differing transmit powers and coverage areas. It would also be of interest to explore how 
various multiple antenna techniques, opportunistic scheduling, and base station cooperation affect 
coverage and rate. 
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Appendix A 



Proof of Theorem [2] 



The proof tracks the proof of Theorem 1 up until step (a) of (PT2i Then, 



ie$\{M 



ft + sR- a 




(39) 
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which admits a much simpler form than (PT21) due to the new assumption that ^ ~ exp(/i). The 
integration limits are still from r to oo and plugging in s = fiTr a now gives 

T 



T+(v/r) a 



vdv 



Employing a change of variables u = (-^x ) results in 

C Ir (fiTr a ) = exp (-7rr 2 Ap(T, a)) , (40) 

where 

poo -, 

Plugging (|40l) into ([111) with v — >■ r 2 gives the desired result. 

Appendix B 
Proof of Theorem [3] 

The ergodic rate of the typical user is r = E[ln(l + sinr)] where the average is taken over both the 
spatial PPP and the fading distribution. Since for a positive random variable X, E[X] = J t>Q F(X > 
t)dt, it follows similar to Theorems CD and [2] that 

r(A, a) = E[ln(l + sinr)] = J e'^E An (\ + j J 2?rArdr 



r>0 ./t>0 



dt27rArdr 



/ e^ Xr2 [ P[/i>r a (a 2 + J r )(e*-l)]dt27rArdr 

Jr>0 Jt>0 

[ e~ wXr2 [ E(exp(- y ur a (a 2 + / r )(e*-l)))dt27rArdr 

Jr>0 Jt>0 

[ e~ wXr2 [ e- <T2 ^ Q(et - 1) /: /r (/ir a (e' - l))dt2vrArdr. 

Jr>0 Jt>0 



From (fT2l) we obtain 



£,>r«( e < - 1)) = exp (-2»A jf (l - ; - ^ ) .*,) 



e < - 1 



exp — 7rAr / — ^dw 

exp H r2(e '- 1)2/ "/,L-TT^ d ^ 



and the proof is complete. 
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Base stations: big dots. Mobile users: little dots. 




10 15 20 



Fig. 1. Poisson distributed base stations and mobiles, with each mobile associated with the nearest BS. The cell boundaries are shown 
and form a Voronoi tessellation. 



Base stations: big dots. Mobiles: little dots. 
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Fig. 2. A regular square lattice model for cellular base stations with one tier of eight interfering base stations. The base stations are 
marked by circles and the active mobile user in the tagged cell by a cross. 



23 



Actual BS locations in a 4G Urban Network 




X coordinate (km) 

Fig. 3. A 40 x 40 km view of a current base station deployment by a major service provider in a relatively flat urban area, with cell 
boundaries corresponding to a Voronoi tessellation. 
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Fig. 4. Probability of coverage comparison between proposed PPP base station model and square grid model with N = 8, 24 and 
a = 4. The no noise approximation is quite accurate, and it can be seen there is only a slightly lower coverage area with 24 interfering 
base stations versus 8. 



24 



Coverage probability for a =2.5, No noise 
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Coverage probability for a =4, No noise 
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Fig. 5. Probability of coverage for a = 2.5 (left) and a = 4 (right), SNR = 10, exponential interference. The proposed model is a 
lower bound and more accurate at lower path loss exponents. 
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Fig. 6. Poisson distributed base stations, no noise, a — 4 with 4 curves corresponding to lognormal shadowing standard deviations of 
0, 3, and 6 dB and Rayleigh fading interference (without shadowing). 
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Coverage probability for a = 3, 6dB LN shadowing, no noise 



-* — Square Grid N=24 
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Fig. 7. Poisson vs. actual vs. grid base stations for a = 3 with LN interference of 6 dB. 



Frequency bands: 4 Frequency bands: 4 




Fig. 8. Left: A spatial reuse of S = 4 is shown for an actual BS development using a greedy frequency allocation. Right: Corresponds 
to a S = 4 reuse for a Poisson base station network with random frequency allocation. The shaded cells use the same frequency. 
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Coverage probability, no noise, a = 4, 8=2 
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Coverage probability, no noise, a = 4, 8=4 
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Fig. 9. Probability of coverage for frequency reuse factors (5 = 2 (left) and (5 = 4 (right). Lower spatial reuse (higher S) leads to better 
outage performance, and we observe that all 3 curves exhibit similar behavior. 




Fig. 10. Average rate of a typical user with SNR = lOdB for both Poisson-distributed and actual base station locations. The average 
rate is maximized when all the cells use the same frequency and hence the complete bandwidth. 



